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Abstract 

Power losses reduction is one of the main targets for any electrical energy dis¬ 
tribution company. In this paper, we face the problem of joint optimization 
of both topology and network parameters in a real smart grid. We consider 
a portion of the Italian electric distribution network managed by the ACEA 
Distribuzione S.p.A. located in Rome. We perform both the power factor 
correction (PFC) for tuning the generators and the distributed feeder recon¬ 
figuration (DFR) to set the state of the breakers. This joint optimization 
problem is faced considering a suitable objective function and by adopting 
genetic algorithms as global optimization strategy. We analyze admissible 
network configurations, showing that some of these violate constraints on 
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current and voltage at branches and nodes. Such violations depend only on 
pure topological properties of the configurations. We perform tests by feed¬ 
ing the simulation environment with real data concerning hourly samples 
of dissipated and generated active and reactive power values of the ACEA 
smart grid. Results show that removing the configurations violating the elec¬ 
trical constraints from the solution space leads to interesting improvements 
in terms of power loss reduction. To conclude, we provide also an electrical 
interpretation of the phenomenon using graph-based pattern analysis tech¬ 
niques. 

Keywords: Distribution feeder reconfiguration; Power factor correction; 
Power losses minimization; Smart grid; Graph-based pattern analysis. 


1. Introduction 


In recent years, the global warming has led people and companies to 
demand for cleaner energy suppliers. Thus, more and more electricity is 
generated from alternative and heterogeneous sources: wind, solar, biofuel, 
and geothermal plants. This phenomenon is called Distributed Generation 
(DG) 12, 14, 231. A Smart Grid (SG) 15] constitutes the improvement of 
a traditional electrical distribution system, which is conceived to overcome 
the problem of the wide diffusion and high penetration of DGs. A SG can be 
seen as an intelligent network able to integrate all users (i.e., producers and 


consumers) with the ultimate purpose of c 


istri 

10 


Duting the electrical power in 
11, 27, 37]. With the advent 


a safe, efficient, and sustainable fashion 
of SGs, the customers of electrical networks become also energy suppliers 
and the load flow in distribution feeders becomes bidirectional. Moreover, a 
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large number of sensors are installed on the network to obtain a complete 


information on the instantaneous status of the in 


that could be exploited for predicting faults 




restructure - information 
38,1. 


The reduction of power losses is one of the main objectives of energy elec¬ 
trical distribution companies. In the literature 2] it is possible to identify 
two mainstream approaches: Power Factor Correction (PFC) and Distributed 
Feeder Reconfiguration (DFR). The PFC tries to reduce the amount of reac¬ 
tive power present in the network in order to (i) minimize the Joule losses, 
(ii) increase the capacity of the network, and (iii) increase the quality of ser¬ 
vice. The DFR, instead, relies on switching a certain number of breakers, 
physically modifying the topological structure of the network and improving 
its operating conditions. In doing so, operational constraints on the net¬ 
work must be satisfied, such as ensuring that no loops are formed and the 
totality of the loads are supplied. Altering the network configuration affects 
the power losses and relieves overload in the network; thus the DFR prob¬ 
lem can be conceptualized as the task of choosing the status of the network 
breakers resulting in the configuration with minimum power losses, yet still 
satisfying the operational constrains. The main drawback of DFR is that it 
results in a complex combinatorial optimization problem, since the status of 
the switches is non-differentiable. This makes the optimization problem re¬ 
lated to DFR very hard to solve. Many researchers have proposed interesting 
solutions in the past. Civanlar et al. Q] propose a heuristic method based 
on a formula that expresses the change in losses between network configu¬ 
rations before and after the reconfiguration process. Moreover in the paper 
the authors suggest a method for filtering configurations that yield lower re- 
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duction in power losses. Other heuristic approaches are presented in Refs. 


22 


24], All such mentioned methods consist in approximate solutions or lo¬ 
cal optimum of the respective optimization problems. To overcome this issue, 


McDermott et ah 


21] use a genetic algorithm providing a good compromise 


between computational burden and quality of the optimization result. More 
recently, novel meta-heuristic methods based on evolutionary optimization 
algorithms are introduced for the same context, showing good experimental 


results 


19 


20j, 


25 


26 


29]. 


One of the main technical difficulties in dealing with the DFR problem 
using evolutionary algorithms is the so-called “radiality constraint”. The 
fulfillment of this constraint makes inappropriate most of the network con¬ 
figurations achievable by switching the available breakers. For this reason, it 
is necessary to conceive a procedure able to select, among all possible con¬ 
figurations, those fulfilling the radiality constraint. A first solution to this 
problem for the SG of ACEA has been proposed by Storti et al. 35|. The 
authors conclude that due to the high complexity of the DFR problem, a 
desirable method should be able to reduce as much as possible the solution 
space, eliminating undesirable switching options a priori. Such configura¬ 
tions are critical with respect to both topological and electrical constraints. 
In 30], the authors propose a heuristic method to compare the admissible 


network configurations in a purely topological manner, facilitating the opti¬ 
mization algorithm in finding the desired solution. All such studies highlight 
the need to identify undesirable network configurations, in order to reduce 
the required convergence time for the faced optimization problem. In fact, 
undesirable configurations, which cause violations of one or more electrical 
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constraints, introduce a significant and unnecessary increase of the compu¬ 
tational demand for the simulation. This is due to the fact that the solver 
tries to perform PFC and DFR using also configurations that are intrinsically 
critical for particular power loads profiles. 

Concerning other works dealing with the Joule losses minimization prob¬ 
lem on real Smart Grids, in Ref. 18] a genetic algorithm with fuzzy logic 
rules is employed to face the optimal power flow problem in the Algerian 
electric network as a Flexible AC Transmission System, while in Ref. [1] is 
proposed a particle swarm optimization method for solving the optimal reac¬ 
tive power dispatch (ORPD) problem on the Algerian electric power system. 


Although not concerning the Jou 
to cite the work of Corsi et al 


e losses minimization problem, it is worth 
4], where the hierarchical voltage control 
system presently applied on the Italian transmission grid is described in de¬ 
tails. Moreover, in Ref. 32] the use of capacitors and static reactive power 
compensators to ensure the voltage stability on the electrical network serving 
the South-West region of France is studied, aiming to develop an advanced 
system for the control of the reactive power compensation. 

In our previous works 


28 




, we faced both the problems of PFC 


and DFR over a portion of the ACEA electrical grid (ACEA is the company 
managing the entire distribution grid of Rome, Italy), using genetic algo¬ 
rithms as optimization strategy. In this paper we elaborate over such studies 
by Erst analyzing how undesirable configurations affect the optimization pro¬ 
cess in terms of both quality of the optimization result and running time of 
the optimization procedure. Then we give an interpretation of the results 
from an electrical point of view, exploiting a graph-based pattern analysis 
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technique. Notably, we study two prototype graphs representing two classes 
of typical network configurations identified in our data. Our work is framed 


in the research area concerning the application of genetic algorithms to 
the joint PFC and DFR problem. As concerns similar works, in Refs. 


ace 


131 , 


[17( | . [9j active power losses minimization is faced by simultaneous capacitor 
placement and feeder reconfiguration by genetic algorithms. With respect 
to our approach, PFC is solved by installing capacitor banks to compensate 
the losses produced by reactive currents, while in the SG of ACEA the PFC 
problem is faced by regulating the power factors of the distributed generators 
in the network. In fact, the telecontrol system of ACEA is capable to modify 
remotely the set points of each renewable energy generator. In Ref. 33] the 


reader can find an interesting review on optimal placement approaches of 
DGs in power systems for optimizing different objective functions (such as 
power losses minimization), including some algorithms based on evolution¬ 
ary computation and swarm intelligence (e.g., genetic algorithms, ant colony 
optimization, particle swarm optimization). Finally, we stress that in the SG 
of ACEA the position of the DGs cannot be relocated, since we are dealing 
with a real network whose physical characteristics cannot be changed. 

The paper is structured as follows. In Section [21 we first describe the 
optimization problem (Section 12.ip and then we introduce the problem of 
admissibility for a network configuration (Section 12.2p . Section [3] provides 
the essential technical details (in Section 13.11) of the electrical distribution 
network under analysis (the ACEA SG) and the related power loss opti¬ 
mization problem (in Section 13. 2p . In Section [4] we analyze the admissible 
configurations and we introduce the concept of constraint compliant conhg- 
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uration. In Section 0 we present and discuss the results of the optimization, 
providing also an electrical interpretation aimed at providing a justification 
for them. Finally, in Section [6] we draw the conclusions pointing at the future 
directions. 

2. Power Loss Minimization Problem 

In this paper, we consider the joint PFC and DFR problem for minimum 
power losses, satisfying constraints on nodes voltage and branches current as 
well as system operating constraints. 

2.1. Optimization Problem 

In this section we formulate the problem of active power losses minimiza¬ 
tion in SGs through PFC and DFR. The problem consists in finding the 
optimal network parameters and the topological configuration that minimize 
the value of the power losses in the network, considering the constraints im¬ 
posed on voltages and currents due to safety or quality of service issues as 
well as physical topological constraints. Consider an admissible set E of the 
network parameters and a suitable cost function J : E —> R that associates 
a real number to each element in E. Formally, the problem consists in mini¬ 
mizing the function J in E. Mathematically, we can express the cost function 
J G [0,1) as follows: 

_ Goss(k) _ Pgen(k) Goad \ 

{} ~PU k) -- ?Wk) ’ u 

where k G E represents an instance of the network parameters, P geri (k) G 
[Coach oo) is the total power generated by all sources, Goad is the total power 
absorbed by the loads, and their difference Goss(k) represents the total power 
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losses in the network. Notice that in the present formulation of the problem, 
-Pioad is independent by the network parameters k because ACEA can provide 
only a description of the loads based on the profiles of active and reactive 
power measured by the meters installed in the network. 

Let us consider a generic SG characterized by n real parameters, m integer 
parameters, and p nominal parameters. We can express the domain of the 
ordinal parameters as: 

A ord = {k ord 6Px Z m : k“' d < k ord < k°^} , (2) 

in which k°t d and k“ d x represent the vectors of the minimum and maximum 
values of the network ordinal parameters, k ord . Concerning the nominal 
parameters, k nom , the domain is a set A nom of all possible admissible elements 
for such parameters: 

d nom = {k nom G Xi x • • • x X p } , (3) 

in which X,; is a generic nominal set with i G {1, The overall domain A 

is defined as A = A ord x A 110111 ; accordingly k = [k ord , k nom ] G A. In order to 
be valid, a solution k must satisfy the constraints on voltages and currents 
defined below: 

B = {k G A : G/ llin < V-(k) < I/. max , i = 1,..., N} 

C = {ke A: |/,(k)| <j = 1, 

where V^(k) is the voltage magnitude of the i-th node for a fixed instance of 
parameters k, N is the total number of nodes, I^ mm and V^ max are the voltage 
limits for the i-th node, while 1 1 3 (k) | represents the current magnitude of 
the j-th branch for a particular instance of parameters k, R the number of 



branches, and finally /" iax the current upper bound for the j'-th branch. The 
definitions given above allow to define the admissible set E as follows: 


E = A n B n C. 


(5) 


Since it is not practically possible to derive expression (pQ) in closed-form 
as a function of k, in the following we will employ a “standard” GA (a well- 
known derivative-free approach) as global optimization algorithm. Moreover, 
since it is also not practically possible to derive closed-forms for Vj(k) and 
Ij{ k) in (ED, we introduce a new function T(k) used in the optimization 
procedure as a measure of the violation of voltage and current constraints. 
Therefore, the constrained optimization problem defined above is actually 
faced by defining the objective function as a convex combination of the fol¬ 
lowing two conflicting terms: 


F(k) = aJ(k) + (1 — a)T(k), 


( 6 ) 


where a G [0,1] is a parameter used to adjust the relative weight of the power 
losses term, J(k), over the constraints violation term, T(k). Thus, our objec¬ 
tive becomes minimizing the function F( k) in the domain A defined through 
(J2D and (ED, instead of J(k) (HD in E ©. Note that (ED is meaningful if and 
only if both J(k) and T(k) vary in the same range, otherwise the optimiza¬ 
tion problem is not well-posed and it is not guaranteed that minimizing F( k) 
in the domain A gives approximately the same result as minimizing J(k) (HD 
in E (ED- 

The function T(k) is defined as follows: 


r (k) = (1 - 0)T 7 (k) + /ftV(k), 


(7) 
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in which fi G [0,1] is a parameter used to adjust the relative weight of the 
violation of current constraints F/(k) with respect to the term related to 
voltages violation Ty(k). In order to make the constraint violation value 
T(k) of the same order as the cost function value J(k), the functions T/(k) 
and IV (k) are defined as follows: 


I>(k)= max {G V (Vfik)/V™™)}, 

*6{0,...,iV} 

r,(k)= m« {G, (/,(k)//"“)}, 


( 8 ) 


where Vfi om indicates the voltage nominal value on the i-th node. The penalty 
functions used in (JH]), that is, Gy(-) and Gi(-) are graphically shown in Figure 
[0 Further details of the optimization procedure can be found in 34], 


2.2. Admissible Network Configurations 

Consider a general SG consisting of several Medium Voltage (MV) feed¬ 
ers, some High Voltage (HV) substations, some DGs, and several loads. In 
order to perform the power loss minimization through DFR, we decided to 
represent the SG as a non-oriented graph Q(N, E), in which N and E are the 
nodes and the edges of the real network, respectively. We introduce Q(N, E), 
the reduced graph of the network, to properly describe all possible system re¬ 
configurations satisfying the topology constraints. The reduced graph of the 
network Q(N, E) does not contain all the information of the original network 
graph Q(N,E), because for our purposes we only need information about 
the connections of different portions of the network and not their detailed 


internal structure. As described in 
through two main steps: 


35]. Q(N,E) is mapped into Q{N,E) 
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nom 




(b) 


Figure 1: Penalty functions: (a) Gy(-), (b) G/(-). 


11 







• The nodes N of the reduced graph are used to model two different types 
of original nodes N. The first one represents nodes at 150kV providing 
the energy balance of the active and reactive power in the SG. In the 
following sections we will refer to it as HV node. The second one can 
represent a single MV real substation connected to loads and DGs, or a 
set of MV substations, powered by only a single HV substation (virtual 
MV). In both cases, we call this kind of nodes as MV nodes. 

• Edges E of the reduced graph are used to model the topology reconfig¬ 
uration. The series of two switches, installed between two consecutive 
MV substations, are mapped into a single edge (virtual breaker) of the 
reduced graph Q. Each edge is associated with a label representing its 
state, i.e., close or open. 


Figure [2] shows an example of the representation of the SG through the 
reduced graph G{N, E). 



Figure 2: Example of a graph Q(N,E) of the simplified network. Yellow circles, gray 
rectangles, dashed blue arrows, and solid red arrows represent MV nodes, HV nodes, open 
status edges, and close status edges, respectively. 

Using the above notation we can introduce the following definitions: 

Def. 1 (Radial Topology Constraint). A network topology satisfies the 
Radial Topology Constraint iff each MV substation is fed by only one HV 
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substation via only one path. 


Def. 2 (Admissible Configuration). A reduced graph Q(N,E) satisfying 
the radial topology constraint is said to be an admissible configuration of the 
network. 


The graph representation Q(N,E) is used to execute an algorithm that 
performs an exhaustive search of all admissible configurations of the net¬ 
work. The details of the automatic procedure are described in Ref. 35]]. The 
output of such a procedure is a list of binary strings (encoding the admis¬ 
sible configurations) having length equals to the number of edges E of the 
reduced network. Each bit represents the state of the corresponding edge 
(virtual breaker). The network topology is specified through a label associ¬ 
ated with the string of bits, spanning the rows of the list of all admissible 
configurations. Because of the nominal nature of the parameter specifying 
the network topology, the objective function becomes non-differentiable and 
the respective DFR optimization problem is very challenging. For this rea¬ 
son, in this paper we use a heuristic method 30] based on the Hamming 
distance between network configurations, to improve the smoothness of the 
objective function with respect to the variation of the nominal parameter 
representing the topological configuration. This technique allows to treat 
during the optimization process the nominal parameter describing t 
figuration as an ordinal parameter. The reader is referred to Ref. 
the details. 


le con- 


301 ] for 
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3. The ACEA SG Pilot Network 


In this work, we consider a portion of the Italian electric distribution 
network managed by ACEA Distribuzione S.pA., located in the west area 
of Rome. The main ACEA goal is the overall improvement of the (i) quality 
of service related to the continuity of electricity distribution, (ii) capacity of 
the network, and (iii) prize of electricity offered to the users. 

3.1. Network Specifications 

The main specifics of the ACEA network are listed below: 

• N .6 Medium Voltage (MV) feeders (n .5 at 20 kV and n.l at 8.4 kV); 

• N .2 High Voltage (HV) Substations; 

• N .76 MV Substations (n.29 at 20 kV and n.47 at 8.4 kV); 

• about 70 km of MV lines (31 km of underground wires and 38 km of 
air lines); 

• about 1200 Low Voltage (LV) user loads; 

• N .5 DGs (n.5 generator sets n.l photovoltaic generator); 

• A. 106 three-phase breakers; 

• N.l TVR (Thyristor Voltage Regulator). 

In each HV substation there is a transformer that converts the voltage from 
150 kV at the primary winding to 20 kV at the secondary winding (HV/MV 
transformer). The cables, the photovoltaic plant, the MV substations, and 
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the TVR are located in the MV portion of the network, whereas the user loads 
and the five generator sets are located in the LV portion of the network. The 
TVR is a series voltage compensation device. It performs a bi-directional 
voltage regulation that maintains the system voltage within specified ranges. 
The bi-directional relation between the input and the output voltage is de¬ 
fined as follows: 

Rout = Rn + iVtapAR, iV ta p G {0, ±1, ±2, ±3}, (9) 

where the values of R n and Rj Ut are expressed in kV and the AR is 0.1 kV. 
The voltage rated value of R n is 8.4 kV. Each MV substation is equipped 
with 2 breakers (switches) that allow to connect the substation with the 
electrical network in different ways. By changing the status of these switches, 
it is possible to modify the topology of the network. The power quality is 
a very important issue in an electrical network, which in turn determines 
the quality of electrical power provided to consumer devices. The correct 
setting of the electrical limits allows operating electrical systems in a safe way 
without significant loss of performance. In order to protect the customers, 
the Authority for Energy and Gas Q has imposed constraints on voltage and 
current to power delivery companies: 

• the instantaneous voltage of all the nodes of the network must be com¬ 
prised in a range of ± 10 % of nominal voltage; 

• the instantaneous current of all the branches of the network must be 
lower than a threshold. 

1 http://www.autorita.energia.it/it/inglese/index.htm 
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These constraints are taken into consideration in the definition of penalty 
functions showed in Figure [0 If the voltage or current, measured at some 
nodes and branches of the network, exceed the admissible range, the value 
of the penalty function increases dramatically. 

3.2. Power Loss Optimization Problem Customization for the ACE A Net¬ 
work 

Using the network specifications described in Section 13.11 we can cus¬ 
tomize the optimization procedure introduced in Section 12.11 In particular, 
we can control the reactive power of the five generator sets through their 
phase parameter, 0. On the other hand, it is not possible to control the reac¬ 
tive power of the photovoltaic generator. Moreover, it is possible to chose the 
iVtap value of the TVR and the topological configuration of the network se¬ 
lecting it from the set of admissible ones, previously determined. The phases 
of the five generator sets, 0i, 0 2 , 03 , 04 , 05 , will be spanned in a real-valued 
range specified by the capability functions of the corresponding generator 
sets. The tap lV ta p of the TVR will be spanned in the discrete (norrned) 
range defined in ([9]). Finally, according to the list of admissible configura¬ 
tions computed using the procedure described in Section 12.21 the network 
topology is specified by the index, N con f, spanning the rows of such a list. 
In particular, in the SG under consideration, the total number of admissible 
configurations is 390. 

Summarizing, a candidate solution vector of the optimization problem 
reads as k = [0 i, 0 2 , 03, 04, 05, iV ta p, iV con f] and technically belongs to the set 
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A defined by the following ranges, 


-0.2 

< 

01 1 02 

< 

0.45 

-0.2 

< 

03 

< 

0.55 

0.0 

< 

04 

< 

0.64 

-0.32 

< 

05 

< 

0.45 

-3 

< 

Atap 

< 

3 

A con f 

G 

{1, ...,390}, 




( 10 ) 


where {1,..., 390} is the set of indexes of all admissible configurations. We 
remind that, due to the heuristic ordering mentioned in Section f7~21 we treat 
the nominal parameter N con f as ordinal during the optimization procedure. 
Moreover, according with the Authority requirements, in order to be valid, 
a candidate solution k must satisfy the constraints on voltages and currents 
defined below: 

B = {k C A : 0.91/“ < VAk) < l.lld nom , i = 1,..., N\ 

( 11 ) 

C={kcA:|/ J (k)|</f ax ,j = l,...,i?}, 

in which N and R represent the total number of nodes and branches of the 
real network, respectively, whereas I4 nom and JJ nax are the nominal value of 
the voltage of the i-th node and the maximum current allowed in the j-th 
wire, respectively. 


4. Analysis of Admissible Configurations 

According to Def. [2] a reduced graph Q (N, E) satisfying the Radial Topol¬ 
ogy Constraint is an admissible configuration for the reference SG. In Section 
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IQ we introduced the problem of minimization of power loss J(k) defined 
in the domain E in terms of the minimization of the convex function F(k) 
in A. It is worth noting that, in order to make the optimization problem 
well posed, both terms of (|6]) must be normalized in the same range (e.g., 
[0,1]). During several preliminary tests, we noticed that for a certain num¬ 
ber of runs of the electrical network simulation, the objective function, F( k), 
assumes values greater than unity, hence violating such a requirement. In 
these situations the GA seems to have difficulties in minimizing the power 
losses, J(k), during the optimization process. According to a preliminary 
analysis of this phenomenon, this behavior seems to be related only to the 
actual value of the topological parameter N coni . In fact, for a certain selec¬ 
tion of the N con { parameter by the GA, the constraint term T(k) seems to 
increase dramatically compared to the J(k) term, causing the power losses 
minimization procedure to fail. To better understand this aspect, in the fol¬ 
lowing section we perform a numerical analysis of admissible configurations 
of the ACEA SG in terms of the violation of electrical constraints and we 
introduce the concept of Constraint Compliant Configurations. 

4-1. Constraint Compliant Configurations 

First we introduce the concept of Constraint Compliant Configurations 
(CCC) through the following definition. 

Def. 3 (Constraint Compliant Configuration). A reduced graph Q(N , E) 
whose nodes satisfy the voltage constraints and edges satisfy the current con¬ 
straint m is said to be a CCC. 


18 


We performed several experiments to numerically show that in the SG 
under analysis there exist admissible configurations that inherently violate 
constraints on voltage at some nodes and/or on current at some branches. 
This means that, independently of the choice of the parameters (pi, 0 2 , 03, 
04, 05, Wap of the network, constraints in (HI]) are violated. Roughly speak¬ 
ing, considering an admissible configuration represented by its reduced graph 
Gi{N , E ), and the corresponding configuration parameter iV C onf = h the value 
of the objective function F( k) associated with it remains of the same order 
of magnitude as the remaining parameters change. To support this assertion, 
we analyze all admissible configurations. In particular, once the configura¬ 
tion parameter iV C onf is set, we perform a random sampling of the subset 
A defined in (fill along all other dimensions. More precisely, we randomly 
choose 2000 points of the subset A \ X. 

The performed analysis highlights two different behaviors. Here, we re¬ 
port the results for two admissible configurations associated with the two 
different behaviors. The results obtained for the two configurations labeled 
32 and 81 are shown in Figure [3] parts (a) and (b), respectively. 

It is worth noting in Figure [3] part (b) that the values of the objective 
function F( k) are distributed in bands due to the influence of the parameter 
Wap- Comparing part (a) and (b) of Figure [3j it can be observed that the 
mean value (shown as a red line in the figure) of the objective function F( k) is 
much higher (lower) than unity in the first (second) configuration. This fact 
shows that, if the parameter Wonf = 32 is selected, the GA will be probably 
unable to make F( k) lower then unity. In fact, the order of magnitude of 
the objective function value depends only on the parameter Wonf and it does 
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Samples Samples 

(a) (b) 

Figure 3: Objective function F(k) computed at 2000 randomly selected points, maintain¬ 
ing constant the parameter (a) N con { = 32, (b) N con f = 81. The mean value is highlighted 
in red, while the standard deviation in dotted black lines. 

not change with the other parameters. Moreover, by computing the ratio 
r) — a/n between the standard deviation and the mean value for the 32-th 
and 81-th configurations we find ij = 0.0104 and rj = 0.0582, respectively. 
This means that in configuration number 32 the variation of the control 
parameters have a lower influence in the minimization of F( k) with respect 
to configuration number 81. 

To better understand this fact, we analyze the different components of 
(|6l) for the 32-th network configuration. The value of the two terms, J(k) 
and T(k), for all the samples and the respective mean values and standard 
deviations are shown in Figure [4] Remember from ([1]) and (J7|) that the term 
J(k) is responsible for the minimization of the power loss, while the term T(k) 
guarantees that a configuration violating the constraints on voltages and/or 
currents is severely penalized by the optimization algorithm itself (insuring 
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a very low fitness value). From Figure [5] it is possible to observe that the 
increase in the objective function is due to the term T(k), meaning that 
for this configuration the constraints (ITT]) are violated, and all parameters 
(except for iV con f) fail to bring the network in a safe condition. We can 
conclude that the 81 -th configuration is CCC, while the 32-th is not. 




Samples Samples 

(a) (b) 

Figure 4: Random sampling of the subset A , maintaining constant the parameter N con f = 
32. (a) Power Loss Term J(k), (b) Constraints Term T(k). 


f.2. Influence of NCCCs on the Objective Function 

As a consequence of what we have empirically shown in the previous 
section, during the optimization procedure candidates that are admissible 
configurations but not CCC will never be solution of the optimization prob¬ 
lem ([6]) . The presence of such configurations in the solution space negatively 
affects the behavior of the optimization procedure, as we will discuss in detail 
in Section [5j 

Here we give a practical demonstration of the negative influence of those 
configurations that are not constraint compliant (shortened as NCCC) in 
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the objective function (]fi|). We know that in the convex combination (15]) . 
the real parameter a is used to adjust the relative weight of the two terms 
of the expression. Moreover, in order for the optimization problem to be 
well posed, both J(k) and T(k) must be normalized in the same range (in 
this case [0,1]). However, in the previous section we have shown that for 
all admissible configurations, the terms J(k) and T(k) are not necessarily 
normalized in [0,1]. In particular, by defining JiuAx(k) as the maximum 
value that the power loss can assume and T M Ax(k) as the maximum value 
for T(k), from (|6|) . we can derive the following expression: 

F(k) = aJ^ik) ■ JMAx(k) + (1 — a)r A r(k) • r M Ax(k), (12) 


in which Jjv(k) and Tjv(k) are the normalized values in [0,1] of the corre¬ 
sponding functions - normalization is implemented by dividing for the max¬ 
imum value. 

In the following, we derive a value for the weighting parameter a, denoted 
as a eq , such that the convex combination in ([T2]) is safely mapped to the unit 
interval. Let us define a eq as the equivalent (i.e., transformed) a coefficient 
that should be used in F(k) if its terms were normalized. Using simple 
algebra (details of the calculations are provided in Appendix A) we can derive 
the analytical expression for a eq : 

«</MAx(k) 


^eq 


(13) 


a^MAx(k) + (1 — a)r MA x(k) 

As an example, for the 32-th configuration considered in the previous sec¬ 
tion, the maximum value of the constraint term is TMAx(k) = 40.38 and for 
the power loss is JMAx(k) = 0.0265. Therefore, assuming that in principle 
we want to give more weight to the minimization of the power loss, i.e., by 
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setting a = 0.9, during the optimization using (TT51) the effective weight is 
set to a eq ~ 0.0059. This means that for all individuals in the GA popula¬ 
tion with the iV^f parameter set to indexes corresponding to NCCCs, the 
algorithm is almost entirely devoted to reducing the value of the constraint 
term instead of the power loss term. However, in the previous sections we 
have empirically demonstrated that, for such configurations, it is practically 
impossible to bring the network in a safe condition, regardless of the set¬ 
ting of all other parameters. This fact has a negative impact on the entire 
optimization process, as we will further discuss in the following section. 

5. Simulation Results 

In this section, we first compare the GA performances when it is used 
to solve the optimization problem presented in Section 12.11 for the electrical 
network presented in Section 13.11 whether in the admissible set A are in¬ 
cluded only CCCs or all admissible configurations (both NCCCs and CCCs). 
Successively, we also give an electrical interpretation of the obtained results. 

5.1. Optimization Results 

Among the 390 admissible configurations described in Section 1X^1 we first 
manually extract all configurations that are CCCs. The result is that 372 
are uniquely classified: 151 are CCCs, while the remaining 221 are NCCCs. 
We want to compare the performances of the GA in solving the optimization 
problem 13.21 when in the solution space are present or not the NCCCs. Let 
us refer to the case with all admissible network configurations as “Experi¬ 
ment 1”, while with “Experiment 2” we refer to the case considering CCCs 
only. In both experiments, we employ the ordering criteria for the network 
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configurations described in Ref. 30] • We use a simulation model of the 


ACEA network 


34] realized using MATLAB and Simulink, together with 


the GA implemented as described in Ref. [7]. To perform all experiments, 
we consider as input of the network model the power profile of distributed 
generators and loads registered at 1:00PM (for one hour) on January 1st. 

In the following, we provide some relevant details about the setting of 
the GA. The number of individuals in the GA population is set to 20; the 
elite individuals are 2, i.e., only 2 individuals in the current generation are 
guaranteed to survive in the next generation; the crossover fraction parameter 
is 0.8; the mutation operator is applied to the remaining individuals with rate 
0.1. Furthermore, the a and f3 coefficients used in expressions (J6]) and ([7]) are 
set to 0.9 and 0.2, respectively. The maximum number of iterations before 
the algorithm halts is 100, but the GA might stop if the relative change of 
the fitness value over 50 iterations is less than or equal to 10~ 9 . We execute 
the GA ten different times with different random initialization seeds; the j- 
th execution considers an equivalent initial population Pj for both series of 
simulations (i.e. only CCCs or all admissible configurations in the solution 
space). However, for simulations that consider CCCs only, the individuals of 
Pj , whose 7-th gene (corresponding to the N coni variable) specifies a NCCC 
are replaced by randomly generated individuals, whose 7-th gene is forced to 
code for a CCC. 

Results of the simulations are shown in Table [H The table reports the 
mean value and standard deviation of the number of generations (^gen) re¬ 
quired for convergence, the fitness value percentage reduction (A F), and the 
reduction of power loss in the network (APi oss ) for both experiment settings. 
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More precisely, the last two indicators compare the fitness value and actual 
power loss at optimal solution with respect to fitness value and power loss of 
the best individual in the initial population. 

Table 1: Mean and standard deviation for number of generations gen) required by the 
GA converge, reduction of the fitness value expressed in percentage (A F), and reduction 
of the power loss (APi oss ) in Experiment 1 (all configurations) and Experiment 2 (CCCs 
only). _ 



Experiment 1 

Experiment 2 

# gen 

73 ± 18.3763 

65 ± 11.6858 

A F [%] 

0.0227 ±0.0001 

0.0205 ±0.0010 

APloss [W] 

1097 ±6.5015 

1138 ± 9.6788 


By analyzing the results in Table [T] we can observe that (i) there is 
no statistically significant difference for the average number of iterations 
required for convergence and (ii) the mean reduction of the achieved fitness 
function value in Experiment 2 is slightly less (but statistically significant) 
than the one observed for Experiment 1. This second fact would led us 
to assume an inferior mean reduction of the power loss in Experiment 2 
as well. However, the power loss reduction achieved in Experiment 2 is 
significantly higher than the one observed for Experiment 1 - please note 
that differences are statistically significant and are evaluated with t-test, 
p < 0.0001. The reason behind this result can be found in the redefinition 
of the solution set considered for the optimization problem in Experiment 
2, composed by CCCs only. In fact, in the set of network configurations 
considered in Experiment 1, many candidate solutions systematically violate 
constraints on voltages and currents causing the problems described in the 
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previous section. Such a behavior is even more accentuated by the fact that 
we have chosen a small value of the coefficient /3, and configurations that 
do not satisfy the constraints cause particularly very strong violations of the 
electrical current constraints. At first, the herein reported overall power loss 
reduction might appear not very significant. However, we remind that all 
tests simulate only one hour of one specific day of the year. Projecting the 
achieved power loss reduction to the entire network over a longer period of 
time might led to important improvements of the operating condition of the 
SG managed by ACEA. 

In conclusion, we demonstrated that removing from the set of admissible 
configurations those that are NCCCs causes a significant reduction of the 
power loss. This fact is particularly important in order to adopt such an 
optimization approach in a real time control system, where desired solutions 
must be determined on a hourly basis, usually relying on limited computa¬ 
tional resources. 

5.2. Electrical Interpretation of the Representative Network Configurations 

As observed in the previous sections, among the admissible configurations 
it is possible to recognize a subset where current or voltage constraints are 
systematically violated (NCCCs), regardless of the DGs setting. 

In order to provide a meaningful electrical interpretation for this fact, a 
single representative network configuration for both sets of CCCs and NCCCs 
is selected. Let S be a set of n graphs (network configurations in our case). 
A natural candidate to represent S is the graph Q* that minimizes the sum 
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of distances (MinSOD) 8], which is determined by the following expression: 

n 

Q* = arg min E d{Gj,Gi). (14) 

Sj€S i= i 


In Eq. [14l is a dissimilarity measure between graphs 16], which 

in our case is implemented as the Hamming distance among the adjacency 
matrix representations of graphs - we remind that we are considering Boolean 
graphs of the same order, i.e., graphs with the same number of nodes entirely 
described by the presence or absence of edges. 



(a) 



(b) 


Figure 5: MinSOD graphs for (a) CCC and (b) NCCC network configuration classes. 

We computed the MinSOD element for the CCC and NCCC classes sep¬ 
arately. The MinSOD graphs are graphically represented in Figure [5] First 
of all, it can be noted that for the representative of NCCCs (J5]part (b)) all 
MV stations are fed by a single HV node, whereas the second available HV 
node is electrically isolated. This configuration results in a very long feeder. 
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As it is well-known, long feeders can suffer of low-voltage problems in the 
nodes far from the feed bar (HV nodes) and can exhibit over-current issues 
close to the feed bar due to the large amount of user load connected to it. 
For this reason, configurations containing very long feeders are very likely to 
be NCCCs, regardless of the DGs setting. Conversely, by analyzing the rep¬ 
resentative graph of CCCs (J5] part (a)) it is possible to note that user loads 
are distributed between the two available HV nodes, resulting in shorter and 
more equilibrated feeders having higher probability to fulfil constraints on 
voltage and current. 

This observation suggests that, for each admissible network configuration 
considered in the optimization problem, the user load should be balanced 
among all the available HV nodes as much as possible, in order to reduce 
the typical length of the feeders and the corresponding amount of required 
electrical current at the feed bar. 


6. Conclusions 


In this paper we 


first described in 
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have presented an improvement of the control system 
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36]. We performed an analysis of admissible net¬ 


work configurations to identify undesirable configurations that may slow 
down the convergence speed of the optimization procedure. In particular, 
we have shown that there exist few admissible configurations for which it 
is not possible to avoid volt age/cur rent constraint violation. We performed 
experiments on real data concerning one hour of power profile of distributed 
generators and loads for the SG located in the west area of Rome, realized 
by the company ACEA Distribuzione. Results showed that, for the network 













under analysis, removing from the solution space of the optimization prob¬ 
lem those configurations that are not constraint compliant leads to overall 
improvements in terms of the power loss reduction. Successively, we used 
graph-based pattern analysis techniques to identify representative networks 
of both the constraint compliant and the undesirable configuration classes. 
We then interpreted the results of such analysis from the electrical point 
of view. In particular, we noted that only those configurations with user 
loads suitably balanced among all the available HV nodes can be actual so¬ 
lutions of the considered optimization problem. Accordingly, configurations 
having very long feeders might be a priori neglected from the search space, 
without causing loss of performance in the optimization of the system. In 
future works, we intend to repeat such analysis for an extended portion of 
the network by considering also different time periods. Moreover, we intend 
to verify if it is possible, using pattern recognition techniques, to predict the 
type of configuration without simulating the entire network, improving thus 
the overall usability of the proposed control system. 

Appendix A. Derivation of Eq. [Ell 

Given the expression ( II21) : 

F(k) = «J N (k) • JMAx(k) + (1 — oOr N (k) • T M Ax(k), 

we want to obtain a normalized expression of the form 

F(k) = a eq J N (k) + (1 — Oeq)rN(k), (A.l) 

in which a eq is the effective parameter of (112]) when both J(k) and T(k) are 
normalized in [0,1] using the respective maximum values. First, we impose 
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that the relative weighting provided by the user-defined a is preserved in the 

analytically calculated a eq . This is done by considering the following ratio: 

qJ N (k) J MAX (k) _ a eq J^(k) 

(1 — a)r N (k)r M Ax(k) (1 — a eq )r N (k) 

Q:T M Ax(k) _ a e q 

(1 — a)r M Ax(k) 1 — a eq 

Accordingly, we can compute a eq by the following manipulations of (1A.2h : 


a 

a 

a 

a 


eq 


= (1 - a 


aTMAx(k) 


eq ' (l-a)r M Ax(k)’ 

«T M Ax(k) _ aJMAx(k) 


eq (1 — a)T MA x(k) eq (1 — a)r M Ax(k) ’ 

(1 — a)r MA x(k) + a ■ JMAx(k) aJ M Ax(k) 


eq 


eq 


(1 — a)r MA x(k) 

«T M Ax(k) 


(1 — a)r MA x(k) 


(1 — a)r MA x(k) + a ■ JMAx(k) 
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